Multiple domain formation induced by modulation instability in two-component 
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The dynamics of multiple domain formation caused by the modulation instability of two- 
component Bose-Einstein condensates in an axially symmetric trap are studied by numerically 
integrating the coupled Gross-Pitaevskii equations. The modulation instability induced by the 
intercomponent mean-field coupling occurs in the out-of-phase fluctuation of the wave function and 
leads to the formation of multiple domains that alternate from one domain to another, where the 
phase of one component jumps across the density dips where the domains of the other exist. This 
behavior is analogous to a soliton train, which explains the origin of the long lifetime of the spin 
domains observed by Miesner et al. [Phys. Rev. Lett. 82, 2228 (1999)]. 
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Spatial pattern formation in systems that are attempt- 
ing to reach equilibrium is one of the most important 
problems in nonlinear and nonequilibrium physics, and 
is encountered in many fields of physical science . The 
modulation instability (MI) is an indispensable mecha- 
nism for understanding pattern formation from a uniform 
system. MI is a general phenomenon in nonlinear wave 
equations, in which weak spatial perturbations with a 
certain range of wave numbers grow exponentially into a 
train of localized waves as a result of an interplay between 
nonlinearity and diffraction effects 0- 

Recently, the MI has been discussed in relation to the 
dynamics of ultra-cold atomic-gas Bose-Einstein conden- 
sates (BECs). This interest was triggered by the exper- 
imental observation of collapsing dynamics P| and the 
generation of matter-wave soliton trains Q in BECs with 
an attractive interaction. The atom-atom interaction in- 
troduces nonlinearity into the system. The key param- 
eter of this interaction is the s-wave scattering length, 
which can be controlled experimentally by the Feshbach 
resonance . Solutions of these dynamics need nonlinear 
analysis beyond the MI; some theoretical and numerical 
studies have reported pattern formation caused by the 
MI in an attractive one-component BEC FR Ula liot . 
Whereas these studies are concerned with the instability 
of self-wave modulation under attractive nonlinearity, if 
there are more than two wave fields, the coupling between 
multicomponent waves, called cross-phase modulation in 
the field of nonlinear optics, can also gives rise to MI 
even if all interactions are repulsive 0, . The present 
study treats the nonlinear dynamics of pattern forma- 
tion caused by cross-phase MI in two-component BECs 
by numerical simulation of the coupled Gross-Pitaevskii 
(GP) equations. 

Pattern formation in a multicomponent BEC was ob- 
served by Miesner et al. Miesner et al. first pre- 
pared a BEC of ^^Na atoms with the hyperfine state 
|F — l,mF = 1) = |1) in a cigar-shaped optical 
trap, and then placed half of the atoms of |1) in the 



\F — l,mp = 0) = |2) state non-adiabatically by us- 
ing rf irradiation. Then, they made use of the quadratic 
Zeeman effect by applying a magnetic field to prevent 
the \F = l,mF = —1) component from appearing. As a 
result, the system could be regarded as binary. Within 
100 msec after the placement, they observed the multi- 
ple domains that alternated from one domain to another. 
These domains had a typical sizes of 40 /im and had a 
very long lifetime of about 20 seconds in the absence of 
a magnetic field gradient. 

Two approaches have been used to gain an understand- 
ing of these long-lived domains. First, mean-field the- 
ory predicts that phase separation occurs in the equi- 
librium state of two-component BECs when the scatter- 
ing lengths satisfy the relation ai2 > t/oio^ |l4| . where 
fli and 02 are the s-wave scattering lengths for same- 
species collisions, and ai2 for different-species ones. Be- 
cause oi = ai 2 = 2.75nm and 02 = 2.65nm for ^^Na 
atoms |l5l Il6| . the condition ai2 > ^^0102 is met and 
the condensates tend to separate spatially. Second, some 
authors pointed out that the multiple domain formation 
is related to MI due to the appearance of imaginary fre- 
quencies in the Bogoliubov modes of two miscible con- 
densates [17I ll^ and analyzed partially the nonlinear 
dynamics after this instability by numerical simulations 
haMml. These finding s are useful for identifying the 
conditions for multiple domains and some of the quali- 
tative features of the dynamics, but there is much more 
to be learned about the nontrivial dynamical features of 
this sytem and the observation in Ref. [l3| . 

This work resolves the above problem by numerically 
solving the coupled GP equations that describe the time 
evolution of two-component BECs at very low tempera- 
tures: 
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with i = 1 or 2, where ipi{r,t) is the condensate wave 
function in the hyperfine state |1) and |2), m the atomic 



2 



mass and fii the chemical potential. The atom-atom in- 
teractions are defined as gi = Anh^ai/m for intracom- 
ponent and gi2 = 4:7rh?ai2/m for intercomponent, corre- 
sponding to self-phase modulation and cross-phase mod- 
ulation in nonlinear optics, respectively. We assume ax- 
ial symmetry (x, y, z) — > (r, z) and use an axisymmet- 
ric harmonic potential Vi = m{u}']_r'^ + uJlz"^) /2. Follow- 
ing the condition of Ref.0], we use the mass of ^^Na 
atoms m — 3.81 x 10~^^ kg, the trapping frequency 
lo±/2ti — 500 Hz and the aspect ratio A = (jj^/cjj^ = 1/70. 
The wave function is normalized by the particle number 
as / dv\^,\^=N^. 

To show that the MI is caused by the cross-phase 
modulation, we first consider the stability of miscible 
two-component BECs in a homogeneous system 0, 0, 
01 , where the equilibrium densities satisfy the relation 
ginio + 5i2^3-z0 = ^-i- When each wave function is as- 
sumed to be written as ipi{Y, t) ~ y^riio+^ie'^'' '"^^^*) with 
the complex- valued amplitude Ai, the linearized equation 
with respect to the fluctuation provides a set of equa- 
tions for the amplitude and gives the dispersion relation 
between the frequency and the wave number k 
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with A = {giTiiQ — (72"-2o)^ + 4(/j2"-io"-20- The condition 
of the MI depends on both the sign and the magnitude of 
the coupling constants. In this work, we focus on the case 
in which the signs of all coupling constants are positive. 
Other cases will be discussed elsewhere. In this case, if 
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+ 5'2"-20 i I (2) FIG. 1: Time development of the condensate densities I'f/'ip 

(solid-curve), |'!/'2p (dashed-curve) and total density nr = 
I^Z-iP + (dotted -curve) along the 2-axis at r = for 

TVi = iV2 = 10^, a-i = ai2 = 2.75nm, 02 = 2.65nm and A = 
1/70. Time of each figure is (a) f = 40 msec, (b) t = 80 msec, 
(c) t — 120 msec and (d) t = 160 msec. The spatial grid is 
[Nr, iVz] = [50, 4400] with the mesh size = = 0.14 fim. 
Figures (e) shows time development of condensate densities 
at r = and z — 0. 



is satisfied, the branch f2?_ becomes negative, which im- 
plies an appearance of purely imaginary frequencies in 
the excitation modes and thus the miscible condensates 
are unstable. Equation Q gives the necessary condition 
312 5i52 for an instability to occur in the range of wave 
number < k < kc = [2m(\/A — 5i?^io ~ 92^20)/^'^]^^'^ ■ 
This instability is of particular importance because it oc- 
curs even for condensates with repulsive interactions, in 
contrast MI in a homogeneous single-component conden- 
sate occurs for only attractive interactions [l]| . 

To see the dynamics caused by the MI, we did nu- 
merical simulations of the coupled GP Eq. The 
numerical method used for solving the time-dependent 
GP equations was the Crank-Nicholson implicit scheme. 
As an initial state, we prepared the stationary solution 
consisting of only the |1) state by solving (— 7i^V^/2m + 
Vi + 5i|V'oP)V'o = MiV'o, where we set iV = / dr\ijjo\'^ = 
2 X 10^. Because the condensate is trapped in a highly 
elongated potential with A — 1/70, the axial direction 
satisfies the Thomas-Fermi limit fii/fiujz — 610 ^ 1. In 
contrast, the radial direction does not satisfy this limit 
because fii/huj± ~ 8.7 0{1). In the experiment, an 



rf pulse was applied instantaneously to transfer half of 
component |1) to component |2). In our simulation, we 
used the same wave function ipo for both initial states |1) 
and 1 2), but their particle numbers are half of the orig- 
inal one, that is, iVi = A^2 = 1 X ID''. The sudden turn 
on of the intercomponent coupling ai2 and the mismatch 
between ai and 02 of the initial configuration result in 
perturbations of the condensates. 

Figure ^ shows the dynamics of multiple domain for- 
mation for oi = ai2 = 2.75nm and 02 = 2.65nm. As 
soon as the simulation began, rapidly oscillating density 
ripples were excited at the central region as shown in Fig. 
n^a), propagating along the z-axis. The density wave is 
out-of-phase between the two components, so that the 
total density tt-t — iV'iP + IV'2p keeps its initial shape. 
This oscillation comes from the radial degree of freedom, 
which is tightly confined so that the mode neither grows 
nor does it form a multi-domain [see also Fig. Qle)]. As 
time passes, because oi > 02, the axial width of |'0i|^ 
expands rather than ip2 as seen in Fig. Qlb). Then, the 
amplitude of the density waves grows first near the edge 
of the condensate [Fig. Qlb) and (c)], eventually devel- 
oping into multiple condensate domains as shown in Fig. 
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^d) . The domains of two components locate alternately 
and the total density does not change even after the do- 
main formation. The domain formation starts at about 
120 msec and the size of one domain is about 20 /ini, in 
good agreement with the experimental observation [l3j. 
Figure ^^e) shows the time evolution of the central den- 
sity \tpi{r = 0,z = 0)p. After the instability has grown, 
the oscillations with a large period and large amplitude 
appear as shown in Fig. QJe), representing the entire 
exchange of the central density of ipi and V'2- 

The growth time of the MI may be determined by 
the fastest growth mode which has the most negative 
value of A straightforward calculation yields fi^^ — 

—h'^kj/Am'^, where kf — kc/V^ represents the wave 

number of the fastest growth mode. To estimate the 

growth time r/, we assume that the density profile of 

the initial condensate |1) has the Thomas-Fermi profile 

noir,z) = IVoP = - (r/Rr? - {z/R ,f)/g^ wi th 

m = (155iA7V/16\/27r)2/5(TOtj2)3/5^ ^ ^2/.ti/mcJi = 

3.9 jim. and Rz = Rr/X = 273 /xm. Because half of the 
condensate |1) is suddenly transferred to |2), we use the 
density no(0,0) with N/2 for mo {i = 1,2) in Eq. 101. 
For iV = 2 X 10^ we obtain r/ = 27r/|f2_/| = 26 msec 
and kf = 0.42 /zm~^. The value of kf corresponds to the 
domain size id = '^T^/kf = 15 yum. Because of Id > Rr 
the instability does not occur in the radial direction. The 
estimated growth time of 26 msec and the domain size 
of 15 fim are in reasonable agreement with the numerical 
results as well as the experimental results. 

It is interesting to point out that the dynamics de- 
scribed here is analogous with the collapse dynamics and 
soliton-train formation in a BEC with an attractive in- 
teraction H, S H 0j ^1 ■ Since the two-component con- 
densates of ^^Na atoms have gi = gi2 = g, Eq. |^ can be 
rewritten as ihdipi/dt = {—h^V^/2m + Vi + gnT)'>pi and 
thd^2/dt = {-n^V^/2m + V2+ griT + (52 - g)\H^)i'2. 
Here, the total density tit = iV'iP + 1^2 P functions as 
the static potential because it hardly changes during the 
time evolution as seen in Fig. Then, an effective at- 
tractive interaction a ppe ars for the condensate |2) be- 
cause 32 - .9 < [H 111. This means that the sudden 
population transfer from |1) to |2) is formally equiva- 
lent to the sudden change of the atomic interaction of |2) 
from positive to negative. Therefore, the spontaneously 
formed domains may have a solitary wave structure like 
that of a bright soliton train 0, S S ^ ■ This is clearly 
seen in Fig. |21 where a spatial distribution of the con- 
densate phase 9i =a,rgipi is plotted for a typical multi- 
domain structure. We can see that the the value of the 
phase of one component is almost flat within each domain 
and jumps across the density dips where the domains of 
the other exist. Also, there is a similarity to the nu- 
merical simulation of bright soliton formation in which 
bright solitons first form at the sharp edge of an initially 
prepared condensate with a rectangular shape 1^, .IQj . 
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FIG. 2: The density profile (solid curves) and phase profile 
(dots) for the multi-domain structure of tpi and i/)2 at t = 225 
msec in the simulation of Fig. 



This is explained by the fact that the wavelength (ampli- 
tude) of self-interference fringes in the initial wave func- 
tion is longer (larger) at the edge of the condensate than 
that at the central part Tol. These fringes can be the 
seed of the modulation, which gets the unstable wave- 
length £d first at the edge. In our simulation, such sharp 
edges are created spontaneously^ as seen in Fig. ^b) and 
(c) , because oi > 02 and the density of ipi is thus pushed 
aside by the repulsion with that of V'2- As a result, the 
density of ipi broadens along the axial direction and de- 
velops a nearly rectangular shape. 

After the multiple domains have formed, the density 
in each domain exhibits the complex dynamics. As ex- 
pected for a bright soliton train in a single-component 
BEC, the dynamics may be controlled by the effective 
interaction between the solitary waves, which depends 
on their phase difference 0, H, 0|. In our simulation, 
the phase difference between domains shown in Fig. |21is 
determined nontrivially, following the nonlinear dynam- 
ics caused by the MI. Then, an exchange of particles be- 
tween two domains of the same component occurs as in 
the Josephson effect, where the domain of the other com- 
ponent has the role of a potential barrier [l^]. The os- 
cillation with a large period and large amplitude shown 
in Fig. ^e) is the cooperative oscillation of the two- 
component soliton trains by this Josephson effect. 

Miesner et al. observed that the multiple domain 
structure persisted for at least 20 seconds 0] , which sug- 
gests that the domain state is relaxed into a metastable 
state by some dissipative process. The reason of the 
metastability was argued to be the large repulsive mean- 
field interaction energy between the two components and 
the large kinetic energy that prevents the wave function 
from speading out in the tight radial direction. To show 
the metastability, we integrate Eq. ^ from the devel- 
oped domain state of Fig. ^d) by introducing a phe- 
nomenological dissipation 7 in the left-hand side of Eq. 
dU as id/dt {i- j)d/dt (7 ~ 0(10-^) ^) under 
the constraint of a fixed particle number. If a metastable 
minimum supports the multiple domains, this structure 
must survive even under dissipation. Figure |2| shows the 
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FIG. 3: The decay of the total energy in the numerical sim- 
ulation of Eq. (0 with the phenomenological dissipation 
7 = 0.025 (solid-line) and 7 = 0.05 (dashed-line) starting 
from the multiple domain state of Fig. Q (d) . Snapshots of 
the condensate densities for 7 = 0.025 are shown at times (a) 
- (d). Inset on right shows the densities of the ground state 
with E = 12.64fiwx 

decay ofthe total energy i? = / dr^^^^ ^Tp*{h^\j'^ 12171+ 

Vi -f^i I '0i P /2) V'i +3 12 1 V'l P I V'2 P and some snapshots of the 
density profile. The simulation shows that the number of 
domains gradually decreases with a monotonia decrease 
of E, and also that the system certainly decays to the 
equilibrium state in which one component occupies the 
middle of the trap and the other sandwiches it [23 • The 
decay time of the domains is very long because the soli- 
tons must move to the condensate edges to vanish. This 
result indicates that the multi-domain state seen in Fig. 
^is not metastable, but instead shows that the solitonic 
character of the multiple domains can produce a long 
lifetime; it is generally known that the decay of the topo- 
logical defects takes a long relaxation time. Although we 
do not know the specific value of 7 13 for Miesner et 
al.^s experiment, the experimental condition had negli- 
gible thermal components and thus the above solitonic 
character can explain the extremely long observed life- 
time of the multi-domain structure. 

In conclusion, we have revealed the nonlinear dynamics 
of multiple domain formation in two-component BECs 
and obtained the excellent agreement with the experi- 
mental results in Ref. [l^ . The underlying mechanism 
of the domain formation is MI induced by cross-phase 
modulation, which occurs in condensates with repulsive 
nonlinearity. This study gives a new interpretation of 
the observed multiple domains, which are two-component 
soliton trains and whose long lifetime is ensured by the 
phase constraint inherent in solitons. Because the pres- 
ence of MI and solitons are interrelated, our results are 



closely related with the soliton excitation in multicom- 
ponent BECs [2^ and may be important for the devel- 
opment of nonlinear atom-optics. 
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